x_max = 2*pi;
eps = 1e-3;
T = 1;
N = 256;
k = 1e-3;
max_iter = round(T/k);
init_u = initial_data1(N,x_max);
ref_u = CrankNicholson_diffusionConvection(init_u,N,x_max,k,4,max_iter);


N1 = 35;
k = 1e-3;
max_iter = round(T/k);
init_u = initial_data1(N1,x_max);
st1=cputime;
u1 = CrankNicholson_diffusionConvection(init_u,N1,x_max,k,2,max_iter);
time1=cputime-st1
err1 = max_norm(ref_u,u1,x_max)


N2 = 37;
init_u = initial_data1(N2,x_max);
st2=cputime;
for i = 1:10
u2 = CrankNicholson_diffusionConvection(init_u,N2,x_max,k,4,max_iter);
end
time2=cputime-st2
err2 = max_norm(ref_u,u2,x_max)

